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Abstract. It is well known that the dominant eigenvalue of a real essentially nonnegative 
matrix is a convex function of its diagonal entries. This convexity is of practical importance 
in population biology, graph theory, demography, analytic hierarchy process and so on. In 
this paper, the concept of essentially nonnegativity is extended from matrices to higher 
order tensors, and the convexity and log convexity of dominant eigenvalues for such a class 
of tensors are established. Particularly, for any nonnegative tensor, the spectral radius 
turns out to be the dominant eigenvalue and hence possesses these convexities. Finally, an 
algorithm is given to calculate the dominant eigenvalue, and numerical results are reported 
to show the effectiveness of the proposed algorithm. 
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1. Introduction 

Tensors are increasingly ubiquitous in various areas of applied, computational, and in- 
dustrial mathematics and have wide applications in data analysis and mining, information 
science, signal/image processing, and computational biology, etc; see the workshop report 
[26] and references therein. A tensor can be regarded as a higher-order generalization of a 
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matrix, which takes the form 

= , Ai^...i^eR, l<ii,...,im<n. 

Such a multi-array A is said to be an m-order n- dimensional square real tensor with 
entries Ai^...i^. In this regard, a vector is a first-order tensor and a matrix is a second-order 
tensor. Tensors of order more than two are called higher-order tensors. 

Analogous with that of matrices, the theory of eigenvalues and eigenvectors is one of the 
fundamental and essential components in tensor analysis. 72 references on eigenvalues of 
tensors can be found in the bibliography [22]. Wide range of practical applications can be 
found the references there. Compared with that of matrices, eigenvalue problems for higher- 
order tensors are nonlinear due to their multilinear structure. Various types of eigenvalues 
are defined for higher-order tensors in the setting of multilinear algebra. For example, the 
eigenvalue, the if-eigenvalue, the i?-eigenvalue, the Z-eigenvalue, the A^-eigenvalue defined 
by Qi for even order symmetric tensors [21], the P eigenvalues for general order symmetric 
tensors, and the mode-? eigenvalues for general square tensors defined by Lim [16], the M- 
eigenvalue for a partially symmetric fourth-order tensor, defined by Qi, Dai and Han [23] , 
the D-eigenvalue for a fourth-order symmetric tensor and a second-order symmetric tensor, 
defined by Qi, Wang and Wu [2l], eigenvalues of general square tensors extended by Qi in 
[22] Chang, Pearson and Zhang in [5] and equivalent eigenvalue pair classes by Cartwright 
and Sturmfels [3]. Here, we are concerned with the one in [5l [22] as reviewed below. 

Definition 1.1 Let C be the complex field. A pair {X,x) G C x (C**\{0}) is called an 
eigenvalue-eigenvector pair of A, if they satisfy: 

Ax""-^ = Ax[™-^], (1) 

where n-dimensional column vectors Ax"^~^ and x^"^~'^'^ are defined as 

Ax := I ^ ] Aii^.-.i^Xi^ ■ ■ ■ Xi^ I and x^ ^ := )i<i<„' 

\l2,...,lm = l / l<i<n 

respectively. 

Nonnegative tensors, arising from multilinear pagerank [16], spectral hypergraph theory 
[HEHH], and higher-order Markov chains [E], etc., form a singularly important class of ten- 
sors and have attracted more and more attention since they share some intrinsic properties 
with those of the nonnegative matrices. One of those properties is the Perron- Frobenius the- 
orem on eigenvalues. In [1], Chang, Pearson, and Zhang generalized the Perron- Frobenius 
theorem for nonnegative matrices to irreducible nonnegative tensors. In [U], Friedland, 
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Gaubert and Han generalized the Perron- Frobenius theorem to weakly irreducible nonnega- 
tive tensors. Further generalization of the Perron- Frobenius theorem to nonnegative tensors 
can be found in [201 123 [2H]- Numerical methods for finding the spectral radius of nonnega- 
tive tensors are subsequently proposed. Ng, Qi, and Zhou [H] provided an iterative method 
to find the largest eigenvalue of an irreducible nonnegative tensor by extending the Collatz 
method [9] for calculating the spectral radius of an irreducible nonnegative matrix. The 
Ng-Qi-Zhou method is efficient but it is not always convergent for irreducible nonnegative 
tensors. Chang, Pearson and Zhang [6] extended the notion of primitive matrices into the 
realm of tensors, and established the convergence of the Ng-Qi-Zhou method for primitive 
tensors. Zhang and Qi ^H] established global linear convergence of the Ng-Qi-Zhou method 
for essentially positive tensors. Liu, Zhou and Ibrahim [17] proposed an always convergent 
algorithm for computing the largest eigenvalue of an irreducible nonnegative tensors. Zhang, 
Qi, and Xu [30j established its explicit linear convergence rate for weakly positive tensors. 

The essentially nonnegative tensor we defined in this paper is ultimately related to the 
nonnegative tensor and includes the latter one as a special case. It is a higher order gen- 
eralization of the so-called essentially nonnegative matrix, whose off-diagonal entries are all 
nonnegative. Such a class of matrices possesses nice properties on eigenvalues. It follows 
from the famous Perron- Frobenius theorem for nonnegative matrices that for any essentially 
nonnegative matrix A, there exists a real eigenvalue with a nonnegative eigenvector, which 
is the largest one among real parts of all other eigenvalues of A. This special eigenvalue, 
termed as r{A), is often called the dominant eigenvalue of A. Moreover, r{A) is known as 
a convex function of the diagonal entries of A. This convexity is a fundamental property 
for essentially nonnegative matrices P [TOl [131 HH] and has numerous applications, not only 
in many branches of mathematics, such as graph theory [25] , differential equations [19] , but 
also in practical fields, e.g., population biology [19], demography [7], and analytic hierarchy 
process as well [12] . A natural question arises: does this convexity maintain for higher-order 
essentially nonnegative tensors? In this paper, we will give an affirmative answer to this 
question. 

Similar to the essentially nonnegative matrix, an essentially nonnegative tensor has a 
real eigenvalue with the property that it is greater than or equal to the real part of every 
eigenvalue of A. We also call it the dominant eigenvalue of A, and denoted by A(^). Par- 
ticularly, if A is nonnegative, we have p{A) = X{A), where p{A) is the spectral radius of A. 
By employing the technique proposed in [19], we manage to obtain that the dominant eigen- 
value is a convex function of the diagonal elements for any essentially nonnegative tensor. 
In addition, it is also a convex function of all elements of a tensor in some special convex 
set of tensors. Furthermore, the log convexity is also exploited for essentially nonnegative 
tensors with whose entries are either identically zero or log convex of some real univariate 
functions. Finally, we propose an algorithm to calculate the dominant eigenvalue, conver- 
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gence of the proposed algorithm is estabhshed and numerical results are reported to show 
the effectiveness of the proposed algorithm. 

This paper is organized as follows. In Section 2, we recall some preliminary results, intro- 
duce the concept of essentially nonnegative tensors, and characterize some basic properties 
of such tensors. In Section 3, we show that the spectral radius of nonnegative tensors is a 
convex function of the diagonal elements, and so is the dominant eigenvalue of essentially 
nonnegative tensors. Section 4 is devoted to the log convexity of the dominant eigenvalue. 
In Section 5, we give an algorithm to calculate the dominant eigenvalue, and some numerical 
results are reported. Some concluding remarks are made in Section 6. 

2. Preliminaries and essentially nonnegative tensors 

We start this section with some fundamental notions and properties on tensors. An m- 
order n-dimensional tensor A is called nonnegative (or, respectively, positive) if Ai^...i^ > 
(or, respectively, Ai^...i^ > 0). The m-order n-dimensional unit tensor, denoted by X, is the 
tensor whose entries are ^j^.. with = 1 if and only if zi = ■ ■ ■ = and otherwise 

zero. The symbol A > B means that ^ — i3 is a nonnegative tensor. A tensor A is called 
reducible, if there exists a nonempty proper index subset J C {1, 2, . . . , n} such that 

= 0, Vzi G /, yi2,...,im^ I- 

Otherwise, we say A is irreducible. We call p{A) the spectral radius of tensor A if 

p{A) = max{|A| : A is an eigenvalue of A}, 

where |A| denotes the modulus of A. An immediate consequence on the spectral radius 
follows directly from Corollary 3 in |21] . 

Lemma 2.1 Let A be an m-order n-dimensional tensor. Suppose that B = a{A + hX), 
where a and b are two real numbers. Then fi is an eigenvalue of B if and only if p. = a{X + b) 
and A is an eigenvalue of A. In this case, they have the same eigenvectors. Moreover, 
p{B)<\a\ {p{A) + \b\). 

Let P := {x e R" : > 0, 1 < z < rz}, and int(P) = {x G : Xj > 0, 1 < z < tz}. 
The Perron- Frobenius theorem for nonnegative tensors is as below, following by [4j Theorem 
L4]. 

Theorem 2.1 If A is an irreducible nonnegative tensor of order m and dimension n, then 
there exist Aq > and Xq G int(F) such that 
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Moreover, if A is an eigenvalue with a nonnegative eigenvector, then A = Aq. If \ is an 
eigenvalue of A, then |A| < Aq. 

The well-known Collatz minimax theorem [9] for irreducible nonnegative matrices has 
been extended to irreducible nonnegative tensors in [U Theorem 4.2]. 

Theorem 2.2 Assume that A is an irreducible nonnegative tensor of order m dimension n. 
Then 

mm max — j — = Aq = max mm — j — , 

xeint(P) Xi>0 zeint(P) Xi>0 xf 

where Aq is the unique positive eigenvalue corresponding to a positive eigenvector. 

For nonnegative tensors, Yang and Yang [28] asserted that the spectral radius is an 
eigenvalue, which is a generalization of the weak Perron- Frobenius theorem for nonnegative 
matrices. We state it [281 Theorem 2.3 and Lemma 5.8] in the following theorem. 

Theorem 2.3 Assume that A is a nonnegative tensor of order m dimension n, then p{A) 
is an eigenvalue of A with a nonzero nonnegative eigenvector. Moreover, for any x G int(P) 
we have 

mm ^ — < p{A) < max ^ — . 

The following inequality and continuity of the spectral radius were given in [28l Lemma 3.5] 
and the proof of [25| Theorem 2.3], respectively. 

Lemma 2.2 Let A be a nonnegative tensor of order m and dimension n, and e > be a 
sufficiently small number. Suppose A< B, then p{A) < p{B). Furthermore, if Ae = A + £ 
where £ denotes the tensor with every entry being e, then 

\iYIip{Ae)= p{A). 

£->-0 

Based on the above results, we can easily get the following lemma. 

Lemma 2.3 Suppose that A is an irreducible nonnegative tensor of order m dimension n 
and that there exists a nonzero vector x G P and a real number /3 such that 

Ax""-' < (3x^"'^'y (2) 

Then /3 > 0, x G int(P), and p{A) < (5. Furthermore, p{A) < j3 unless equality holds in 



5 



Proof. Assume on the contrary that for x G int(P) there exists a nonempty proper index 
subset Jc{l,2,...,n} such that Xj = for « G / and > for i ^ I. It follows from ([2]) 
that 

= 0, Vzi G /, VZ2, . . . ^ /. 

A contradiction to the irreducibility of A comes, which henceforth implies that x G int(P). 
Together with Lemma 2.2 in [TS], Ax'^^^ G int(P) is established. It further deduces that 
/3 > 0, and then the last statement holds from Lemma 5.9 in [28]. This completes the proof. ■ 

A simple but useful result follows immediately from Lemmas 12.21 and 12.31 

Lemma 2.4 Let A and E he irreducible nonnegative tensors of order m dimension n. If 
A<B and A^ B, then p{A) < p{B). 

Proof. By Lemma I2l2l p{A) < p{B). Since B is irreducible, Theorem 12 . 1 1 implies that there 
exists X G int(P) such that 

Ax""'^ < Bx""'^ = p{B)x^"'-^\ (3) 

Since x G int(P) and A B, equality cannot hold in ([3]). The desired strict inequality 
p{A) < p{B) holds from Lemma 12.31 ■ 

The remaining of this section is devoted to the essentially nonnegative tensor, with the 
introduction of its definition and some basic properties. 

Definition 2.1 Let A he an m-order and n-dimensional tensor. A is said to he essentially 
nonnegative if all its off-diagonal entries are nonnegative. 

Theorem 2.4 Let A he an m-order and n-dimensional essentially nonnegative tensor. Then 
there exists a > such that aX + A is nonnegative. Moreover, A has a real eigenvalue 
X{A) with corresponding eigenvector in P and X{A) > ReA for every eigenvalue A of A. 
Furthermore, 

\[A) = p{al -^A) - a. 

Proof. Take 

a = max \Ai A + 1. 

l<i<n 

Clearly, a > and aX + ^ is nonnegative. By Lemma 12.11 and Theorem 12.31 we have 

p(aX + ^) = a + Ai, (4) 
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where Ai is an eigenvalue of A with corresponding eigenvector in P. Thus, (jl]) imphes 
Ai G R. Let A(^) = Ai, It follows from Lemma [2. II that. 



X{A) + a = max{|a + A| : A is an eigenvalue of A} 
> la + A| > a + ReA. 

The desired result arrives. ■ 

We call such an eigenvalue in the above theorem the dominant eigenvalue of A. Through- 
out this paper, p{A) and A(^) will denote the spectral radius and dominant eigenvalue re- 
spectively of a tensor A. In the next section, we will show that both p{A) and A(^) are 
convex functions of the diagonal elements of A. 



3. Convexity of the spectral radius and the dominant eigenvalue 

Based on Theorems 12.11 and 12 . 3[ we proceed with the convexity of the dominant eigenvalue 
of essentially nonnegative tensors in this section. It can be verified that the diagonal entries 
have nothing to do with the irreducibility of a tensor. Specifically, let A be an essentially 
nonnegative tensor of order m and dimension n, define a nonnegative tensor B by Bi^,,,i^ = 
a ii = ■ ■ ■ = im and the others are Aj^ . j^. Then A is irreducible if and only if B is. 
Equivalently, A is irreducible if and only if ^ + aX is, whenever it is nonnegative. Thus, by 
Lemma 12.21 and Theorem 12. 4[ it is sufficient to consider the class of irreducible nonnegative 
tensors. 

Theorem 3.1 If A is a given irreducible nonnegative tensor of order m and dimension n, 
andV is allowed to vary in the class of nonnegative diagonal tensors, then the spectral radius 
p{A + V) is a convex function of the diagonal entries ofV. That is, for nonnegative diagonal 
tensors C and T> we have 

p{A + tC + {l-t)V) <tp{A + C) + {l-t)p{A + V), VtG[0,l]. (5) 

Moreover, equality holds in ^ for some t G (0, 1) if and only if V — C is a scalar multiple 
of the unite tensor X. 

Proof. Since both A + C and A + V are irreducible nonnegative tensors, by Theorem 12.11 
and Theorem 12.31 we have p{A + C) > 0, p{A + V)>0, and there exist x,y E int(P) such 
that 

{A + C)x'^'^ = piA + C)a;['"-^1 , {A + V)y"'~' = p{A + V)y^'^-^^ . 
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That is, for i = 1, 2, . . . , n we have 

n 

p{A + C) = Ci,„i+ 

n 



'■112. ..Im 
22...lm = l 



and hence p{A + C) — Cj..i > and pi^A + V) — Dj . j > 0. The inequahty between geometric 
and arithmetic means yields 



i-t 

>...im = l / \t2..im=l ^* / 

+ (l-t)(p(^ + P)- A....). (6) 
Therefore, Holder's inequality and Theorem 12.21 give from ([6]) 

piA + tC + il~t)V) < max \tCi...i + il-t)Di,„,+ V Aii,,„i 

l<i<n ' 

k i2...im=l 

< tp(^ + C) + (l-t)p(^ + P), 

where Zi = xly]"^ for i = 1, . . . , n. This shows ([5]) holds. 

The inequality between geometric and arithmetic means implies that equality in (|5]) holds 
for t e (0,1) if andonly if p(^+C)-Ci...i = p{A+V)-Di,„i for i = 1, . . . , n, i.e., V-C = 7X 
where 7 = p{A + V)— p{A + C). This completes the proof. ■ 

The convexity involved in Theorem 13.11 can be extended to the case of essentially non- 
negative tensors as follows. 

Corollary 3.1 If A is a given irreducible essentially nonnegative tensor of order m dimen- 
sion n and V is allowed to vary in the class of diagonal tensors, then the dominant eigenvalue 
X{A + V) is a convex function of the diagonal entries ofD. That is, for diagonal tensors C 
and T) we have 

X{A + tC + {l-t)V)<t\{A + C) + {l-t)\{A + V), VtG [0,1]. (7) 

Moreover, equality holds in ^ for some t G (0, 1) if and only if V — C is a scalar multiple 
of the unite tensor X. 



Proof. Take 



a = 1 + max{|v4j...j| + |Cj...i| + |A..i|}- 

l<i<n 
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Then dZ + A + C and aX + A + V are all irreducible nonnegative tensors. By Theorem 12.41 
and Theorem 13.11 we have for < t < 1 



Invoking the continuity presented in Lemma 12.21 it is easy to see that Theorem 13.11 and 
Corollary 13.11 hold even when A is reducible. Moreover, Theorem 13.11 and Corollary 13.11 
give necessary and sufficient conditions for the strict convexity. It is worth pointing out 
that the convexity of the dominant eigenvalue only works on the diagonal elements other 
than on all elements of the essentially nonnegative tensor, unless for some special cases. By 
collecting all symmetric essentially nonnegative tensors of order m and dimension n, we can 
get a closed convex cone, says S{m,n). The dominant eigenvalue of any tensor in S{m,n) 
remains convex of all elements of the corresponding tensor in the domain S{m,n), as the 
following proposition shows. 

Proposition 3.1 For any A, B E S{m, n), and any t G [0, 1], we have 



Proof. For any A, B E S{m,n), there exists an integer k > such that A + kX and 
B + kX are nonnegative and symmetric and hence for any of their convex combinations. The 
Perron- Frobenius theorem then ensures that p{A + kX), p{B + kX) and p{tA+ (1 — t)B + kX) 
(t G [0, 1]) all act as eigenvalues of the corresponding nonnegative symmetric tensor. By the 
variational approach, it follows that 



\{A + tC + {l-t)V) + a 



= p{aX + A + tC + {1 - t)V) 

< tp{aX + A + C) + {l-t)p{aX + A + V) 

= tX{A + C) + {l-t)X{A + V) + a, 



which yields This completes the proof. 



X{tA + (1 - t)B) < t\{A) + (1 - t)\{B). 



p{tA + (1 - t)B + kX) 




= tp{A + kX) + {l-t)p{B + kX). 



Combining with the fact that p{A + kX) = \{A) + k, the desired convexity follows. 
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4. Log convexity of the spectral radius and the dominant eigenvalue 



If a function /(x) is positive on its domain and log/(2;) is convex, then f{x) is called log 
convex. It is known [15] that the sum or product of log convex functions is also log convex. 
In this section we extend Kingman's theorem [TJ] for matrices to tensors. Our motivation 
for the following proof comes from [T9] . 

Theorem 4.1 Fort G [0,1] assume that J^{t) = (Fj^,.,j^(t)) is an m- order n- dimensional 
irreducible nonnegative tensor, and suppose that for 1 < ^i, . . . , im < '^j Fh...im{'t) is either 
identically zero or positive and a log convex function of t. Then p{J^(t)) is a log convex 
function of t for t G [0,1]. That is, if J^{0) = A, J^{1) = B, and a nonnegative tensor 
Q{t) = {A^;'..^Bl..J,then 

pim) < pm)) < PiAy-'piBY. (s) 

Moreover, the first equality occurs in ^ for some t with t G (0, 1) if and only if 

m = G{t), 

and the second equality occurs in ^ for some t with t G (0, 1) if and only if there exists a 
constant a > and a positive diagonal matrix D = diag{di, . . . , dn) such that 

m— 1 

B = aA- ■ with = aA,,,,„,y-^'^-''>d,, ■■■d,^. 

Proof. Clearly, ^(0) = J^(0) = A and Q{1) = J^{1) = B. The log convexity assumption on 
Fi^.-.i^if) implies that, for t G [0, 1], 

which, together with Lemma [2. 2 [ implies 

p{:F{t)) < p{g{t)). (9) 

Since J^{t) is irreducible, if equality holds in IH] for some to with < to < Lemma 12.41 
implies that J-'(to) = Qito). 

Since J^(0) and J^(l) are irreducible nonnegative. Theorem 12.11 shows that there exist 
x,y E int(P) such that 

Ax""-' = p(^)a;["^-^] , By""-' = p{B)y^'^-^^ . 

For a fixed t G (0, 1), define z = x^~^y^, i.e., Zi = x\~^y\ for \ < i < n. Then, the zth 
component of Q{t)z^~^ satisfies 

n 

«2---jm = l 
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Hence, Holder's inequality gives 

i2---im = l / \i2...im = l / 

= p{AY''piByzr-'. (10) 

It follows from Lemma 12.31 and fllUp that 

piGit)) < p{Af-'p{Br. 
Furthermore, equality holds in [10] for some t G (0, 1) if and only if, for 1 < z < n. 

Summing (|TT|) over 12 ■ ■ - im yields 

p{B)yr' = o.p{A)xT-\ (12) 

Take 



^ p{Ay m 



Then, combining f|TT]) and f fT2|) we obtain 



R. . . — - • ■ . . 

i.e., 

m— 1 

This completes the proof. ■ 

By Theorems 12.31 and 12.41 the above theorem also holds for the dominant eigenvalue of 
J-'(t), when J^{t) is essentially nonnegative with t G [0, 1]. 

5. An algorithm for calculating the dominant eigenvalue 

Let A be an essentially nonnegative tensor of order m and dimension n. In this section 
we propose an algorithm to calculate the dominant eigenvalue of an essentially nonnega- 
tive tensor. This algorithm is a modification of the Ng-Qi-Zhou algorithm given in [T8] . 
By Lemma 12.21 and Theorem 12. 4[ we modify the Ng-Qi-Zhou algorithm such that for any 
essentially nonnegative tensor, the sequence generated by the modified algorithm always 
converges to its dominant eigenvalue. 
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Define two functions from int(P) to P: 

p (x) := mm ^ — , G x := max ^ — , (Id) 

where W is an irreducible nonnegative tensor. The details of the modified algorithm are 
given as follows. 

Algorithm 5.1: 

Step 0. Given a sufficiently small number e > 0. Let 

W = A + aI + £, (14) 

where 

a = max \Ai A + 1, 

l<i<n 

and £ is the tensor with every entry being e. Choose any x^^^ G int(P). Set 
y(0) = n;(a;(0))'"-'andfc:=0. 



Step 1. Compute 



I- m — 1 - 



According to f[T51) . compute F{x'^^~^^^) and G{x'^''~^^'^). 

Step 2. If G'(x*^'^+^^) — < stop. Output e-approximation of the dominant 

eigenvalue of A: 

^(k+i) ^ 1 ^G{x^k+i)) _^ ^(^(/c+i))) _ (15) 

and the corresponding eigenvector x^''~^^\ Otherwise, set k := k + I and go 
to Step 1. 

Clearly, the tensor W defined by (IMI) is positive and hence it is primitive. By Theorems 
12.11 and 12.21 Algorithm 5.1 is well-defined. As an immediate consequence of Lemma 12.21 
Theorem 12. 4[ and Theorem 3.3 in [20], we have the following convergence theorem. 



Theorem 5.1 Let A be an essentially nonnegative tensor of order m and dimensional n, 
and let W he defined by [T^ where e is a sufficiently small number. Then the sequences 
{F{x^''^)} and {G{x^''^)}, generated by Algorithm 5.1, converge to where Ag is the unique 
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positive eigenvalue ofW. Moreover, the sequence {x^^-*} converges to x* and x* is a positive 
eigenvector of W corresponding to the largest eigenvalue . Furthermore, 

limAe = A*, limx* = x* , 

where A* is the spectral radius of A + aX and x* is the corresponding eigenvector. In 
particular, the dominant eigenvalue of A is X{A) = X* — a and x* is also the eigenvector 
corresponding to A(^). 

Proof. It follows from f lT4|) that W is positive, and hence it is irreducible. Therefore, for 
any nonzero x G P, we have yVx"^~^ G int(P), which shows that the tensor W is primitive. 
Hence, by Theorem 3.3 in pU] . 

lim F{x^''^) = lim {G(x('=)) = A„ lim x'-''^ = x*. 

Therefore, A^ — a is an ^-approximation of the dominant eigenvalue of A from Theorem 12.41 
Furthermore, It follows from Lemma 12.21 that 

limAe = A*, limx* = x* . 

It is easy to see that A* — a is the dominant eigenvalue of A with corresponding eigenvector 

X*. ■ 

The above theorem shows that the convergence of Algorithm 5.1 is established for any 
essentially nonnegative tensor without the irreducible and primitive assumption. In order 
to show the effectiveness of Algorithm 5.1, we used Matlab 7.4 to test it on the following 
three examples. 

Example 5.1 Consider the 3-order ?>- dimensional essentially nonnegative tensor 



where 



A 


= [^(1,:, 


:),A(2,: 


,0,^(3,:,:)], 






/ -1.51 


8.35 1.03 \ 




,1) = 


4.04 


3.72 1.45 






^ 6.71 


6.43 1.35 y 






/ 9.02 


0.78 6.89 \ 




,2) = 


9.71 


-5.32 1.85 






^ 2.09 


4.17 2.98 y 






/ 9.55 


1.57 6.91 \ 




,3) = 


5.63 


5.55 1.43 






^ 5.76 


8.29 -0.15 y 



13 



Example 5.2 Let a 3-order 3- dimensional tensor A be defined by Ai^^ = A233 = A^u = 
A322 = 1; ^111 = ^222 = — 1 md zero otherwise. 

Example 5.3 Let a 3-order A- dimensional tensor A be defined by Am = A222 = ^333 = 
A444 = -1, Au2 = Am = Ai2\ = Aisi = A212 = ^332 = ^443 = 1, and zero otherwise. 

Clearly, the essentially nonnegative tensors defined as Examples 15. II and 15.21 respectively are 
irreducible. While, the essentially nonnegative tensor defined as Example 15.31 is reducible. 
We take e = 10~^ and terminate our iteration when one of the conditions — F(a;*^'^^) < 

10~^ and k > 100 is satisfied. Algorithm 5.1 produces the dominant eigenvalue A(^) = 
36.2757 with eigenvector x* = (1.0000; 0.8351; 0.9415) for Example 15.11 the dominant eigen- 
value A(^) = 1 with eigenvector x* = (0.5000; 0.5000; 1.000) for Example 15. 2^ and the domi- 
nant eigenvalue A(^) = 0.8225 with eigenvector x* = (1.0000; 0.7408; 0.9714; 0.5330) for Ex- 
ample [5]3l The details of numerical results are reported in Tables 1 and 2. We list the output 
details at each iteration for Example 15. II in Table 1. We also report the number of iterations 
{No. Iter), the elapsed CPU time {CPU(sec)), the lower bound A^''^ = F(x(^)) - a and the 
upper bound A^''^ = G{x^''^) - a for A; > 1, the error A^'^) = \\A{x^^^)"'-'^ - A('=)(xW)[™-il ||oo, 
and the approximation A'''^^ defined by f|T5l) of the dominant eigenvalue in Tables 1 and 2. 

From Tables 1 and 2, we see that the sequence generated by Algorithm 5.1 converges to 
the dominant eigenvalue of the essentially nonnegative tensor without irreducibility. Algo- 
rithm 5.1 is promising for calculating the dominant eigenvalues of the test three examples. 



Table 1: Detailed output of Algorithm 5.1 for Example 15.11 



k 






A(^) 






1 


35.9969 


36.5635 


36.2802 


0.5666 


0.2833 


2 


36.2554 


36.3030 


36.2792 


0.0476 


0.0211 


3 


36.2747 


36.2776 


36.2762 


0.0030 


0.0015 


4 


36.2757 


36.2758 


36.2757 


9.1725e-005 


4.5870e-005 


5 


36.2757 


36.2757 


36.2757 


6.7568e-006 


2.9868e-006 


6 


36.2757 


36.2757 


36.2757 


4.6425e-007 


2.2441e-007 


7 


36.2757 


36.2757 


36.2757 


1.9041e-008 


1.4348e-008 


8 


36.2757 


36.2757 


36.2757 


8.8998e-010 


8.1036e-009 
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Table 2: Output of Algorithm 5.1 for Examples 15. H I5.2i and 15.31 



Example 


No.Iter 


CPU(sec) 










A(fc) 




8 


0.013 


36.2757 


36.2757 


36.2757 


8.8998e-010 


8.1036e-009 




31 


0.035 


1.0000 


1.0000 


1.0000 


9.6831e-010 


4.1210e-009 




37 


0.078 


0.8225 


0.8225 


0.8225 


7.3324e-010 


1.0635e-008 



6. Conclusions 

In this paper, we have introduced the concepts of essentially nonnegative tensors, which 
is closely related to nonnegative tensors. The main contribution is the convexity and log con- 
vexity of the dominant eigenvalue of an essentially nonnegative tensor, and hence the same 
for the spectral radius of a nonnegative tensor. By modifying the Ng-Qi-Zhou algorithm 
|18] . we have proposed an algorithm (Algorithm 5.1) for calculating the dominant eigen- 
value. Its convergence can be established for any essentially nonnegative tensor without the 
assumptions of irreducibility and primitiveness. Numerical results indicate that Algorithm 
5.1 is promising. 
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